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Abstract 



We study the Causal Dynamical Triangulation (CDT) with extended interactions in 1+1 dimensions 
applying the method in the non-critical string field theory (SFT) constructed by Ishibashi and Kawai. 
For this model, we solve Schwinger-Dyson's equation (SDE) for disk amplitude perturbatively, and find a 
matrix model in the continuum limit reproducing the SDE in the non-critical SFT approach as the loop 
equation. 
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1 Introduction 



As a consequence of general relativity, "un-countable" lots of physics about Universe have been uncovered. 
Now, we have to go beyond the theory and into the quantum realm, i.e. quantum gravity. However, it 
has been known that there is a difficulty in the case that we extend general relativity to quantum gravity. 
Namely, general relativity is not renormalizable at least perturbatively. 

As a candidate to overcome such a serious problem, a kind of non-perturbative method has been 
proposed, which is called Euclidean Dynamical Triangulation (EDT). In EDT, discretizing Euclidean 
space-time by simplices having the lattice spacing a as each side length, we can carry out the Euclidean 
gravitational path-integral non-perturbatively. An important point here is that the lattice spacing a 
is about the inverse energy cut-off A grav . Unfortunately, in EDT, no reasonable classical space-time 
has been found in 4 dimensions, and what has been found are only the skinny polymer-like geometry 
or the dense crumpled geometry, which has been calculated with the help of Monte-Carlo simulations. 
This is because the geometries based on EDT is too "wild" to handle. Alternatively speaking, infinite 
numbers of baby universes are produced in this approach. However, EDT had a great deal of success in 
2 dimensions. In the suitable continuum limit, physical quantities such as several critical exponents and 
correlation functions in EDT realize those of the so-called quantum Liouville theory. Furthermore, a dual 
expression of EDT has been found, and it is called the matrix model. Utilizing the powerfulness of the 
matrix model, conformal matters realized in the so-called minimal model have been successfully included 
in EDT approach. 

In this line of thought, a kind of breakthrough has been casted out, which is known as Causal 
Dynamical Triangulation (CDT) [1 . In CDT approach, the path-integral of dynamically triangulated 
geometries can be done non-perturbatively under the two new additional restrictions. First, one gives the 
Lorentzian signature to simplices. Second, the time-foliation structure is imposed. In this approach, our 
de-Sitter universe can be "realized" in 3 + 1 dimensions [2]. Furthermore, in 1 + 1 dimensions physical 
quantities, say disk amplitude, can be solved analytically pQ. An outstanding feature of the pure CDT is 
that no baby universe is allowed according to the non-anomalous scaling dimension of time. Related to 
the fact above, for instance the Hausdorff dimension du in the (l+l)-dimensional setup is not anomalous, 
dn = 2, compared to that in EDT, dn = 4. This is one of attractive traits of CDT. 

The CDT approach really restricts the configurations of geometries to the causal ones a priori, but in 
fact we do not understand whether or not we should exclude the contributions from the baby universes and 
furthermore from the geometries with different space-time topologies. Focusing on the (l+l)-dimcnsional 
case, CDT has been extended to the one including topology changing processes within the criterion that 
the scaling behavior does not change, i.e. the causal geometries are still dominant, via the non-critical 
String Field Theory (SFT) In addition, the matrix model expression for the non-critical SFT based 
on CDT has been found [4]. Such extended models including baby universes and topology changes are 
called Generalized CDT's (GCDT). As for the matter-coupled CDT's, there is not any analytical tool 
to calculate even in 1+1 dimensions. From the lessons based on the subsequent works in CDT, it can 
be said that the dominance of causal geometries, characterized by the fact that space and time have 
the same scaling dimension, prevents the stampede of geometries. Alternatively speaking, the causality 
makes geometries obedient to handle. If the CDT approach is on a correct direction as quantum gravity, 
this may be a pretty nice property. 

In this paper, to read off some hidden traits of CDT, we quest for possibilities to extend the GCDT 
approach without changing the scaling dimensions of space and time in 1+1 dimensions. We actually 
extend GCDT applying the method in the non-critical SFT techniques in [5] and [5J. We solve the 
Schwinger-Dyson's equation (SDE) for disk amplitude in our model by the perturbation w.r.t. the string 
coupling constant. Moreover, we define the corresponding matrix model in the continuum limit. In 
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Section 2, we review known facts for GCDT. In Section 3, our extended model is explained in detail. 
Both sections are almost separated by the two different subsections, Non-critical SFT Approach and 
Matrix Model Approach. At the end of Section 3, as a consistency check, we also consider the inclusive 
process, which turns out to reproduce our differential equation for disk amplitude. In Section 4, we 
discuss our model from the two different field theories, the non-critical SFT and the matrix model. 



2 Generalized CDT 

2.1 Non-critical SFT Approach 

We shall review the non-critical SFT of the original GCDT formulated in [3] . This model really reproduces 
the disk amplitude derived in the continuum limit of the strictly causal CDT in the case that the string 
coupling constant is zero. In this model, closed strings with length L are created and annihilated from 
the vacuum, |0) ((0|) by the operators, ij. >'(L) and i)(L), respectively: 

(0|^(i) =^(L)|0)=0. (2.1) 

These creation and annihilation operators obey the following commutation relations: 

[i>(L),^(L')]=S(L-L'), (2.2) 

and the others are zero. The world-sheet which closed strings sweep out can be seen as the whole 
space-time itself. Corresponding Hamiltonian can be written as: 



H = dL^(L)U (L,A)i,(L) + G s dL x / dL 2 ^ '(L x )^ {L 2 )^{L X + L 2 ){L X + L 



-2) 
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+ aG s dL x dL 2 ^(L x + L 2 )ip(L 2 )i/j(L x )L 2 L x - dL5(L)^{L), (2.3) 
Jo Jo Jo 

where 

H (L,A) = -Ldl + AL. (2.4) 

G s and A are the string coupling constant and the cosmological constant, respectively. The parameter 



t 



Propagator Branching interaction Merging interaction Tadpole 

Figure 1: Terms in the Hamiltonian 



a in (2.3) was introduced to count the numbers of genus in amplitudes. In the following discussion we 



shall take a — 0, which suppresses the creation of handles. The Hamiltonian above has been determined 
under the following scaling dimensions: 



IS] 



[^(L)] 



(2.5) 
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where a is the scaling dimension of space, or alternatively speaking the lattice spacing, and [S] is the scal- 
ing dimension of time. A crucial difference between the Hamiltonian of the non-critical SFT constructed 
by Ishibashi and Kawai [5j and that of GCDT is the existence of the propagator term, J dLip^ (L)Hoip(L). 
In GCDT the propagator term actually exists but IK's theory does not. This difference comes from the 
fact that both theories have quite different definitions of "time" . 

The authors in [5] derived Schwiner-Dyson's equation (SDE) for the Laplace-transformed disk ampli- 
tude, W A {Z) =/ oo dLe- LZ (0|e" SH «7/>t(L)|0)| s ^oo, in GCDT as[] 

d z [{A~Z 2 )W A (Z) + G s Wl(Z)]+l = 0. (2.6) 

The solution of the above SDE was also derived by a perturbative expansion w.r.t. the string coupling 
constant in [3J: 

Wa(Z) = l -^-G s Z + 3 ^ +0(G 2 S ). (2.7) 

The first term in the solution above is equivalent to the strictly causal solution lj. In this formalism, 
the contributions from baby universes are weighted by the string coupling constant G s . 

2.2 Matrix Model Approach 

The hermitian N x N matrix model reproducing the SDE of GCDT was introduced. We start with the 
following matrix integral [J]: 

" ^e~f= v ^\ (2.8) 



where 

F(0) = -#+^ 2 -^ 3 , (2.9) 

and (f>, g and g s are a N x N hermitian matrix, the 'tHooft coupling constant and the string coupling 
constant, respectively. Then, by introducing the infinitesimal lattice spacing a, we can expand the 
coupling constants and the matrix w.r.t. a: 

g s = ~a 3 G s , ^ = I-a<S> + 0(a 2 ), g = ±(l - l -a 2 K + 0(a 4 )) , (2.10) 

where / is the unit N x N matrix, and G s , and A are the corresponding renormalized values. Substi- 
tuting the fine-tuned values above into the potential j-V(4>), we find 

—trV(6) = 7=-tr( -$ 3 - A$ )+ (terms independent of $)+0(a). (2.11) 
9s Gs V 3 / 

Here we define the partition function in the continuum limit as: 

N fl 



Z = d$ exp 



-— tr( -$ 3 - A$ 
G s 1 3 



(2.12) 



In the large- A'' limit, the saddle-point equation become^] 

d z [(A - Z 2 )W A (Z) + G S W A (Z) 2 ] +1 = 0, (2.13) 

1 The authors derived the more general result with arbitrary a, but here we restricted our situation to that with a = 0. 
2 In [4], the authors derived the general saddle-point equation beyond the large- N limit. The general saddle-point equation 
really coincides with the SDE with arbitrary a by the treatment, a = 1/N 2 . 
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where Wa(Z) is the resolvent for the matrix <J>. We notice that the saddle-point equation coincides with 
the SDE of GCDT. 



3 Generalized CDT with Extended Interactions 
3.1 Non-critical SFT Approach 

Applying the method in [S] , we shall construct the non-critical SFT Hamiltonian of GCDT with extended 
interactions. 



The propagator term in (2.3), J dLi/j^ (L)Hq(L, A)ifj(L), induces the strictly causal geometry. To make 



this propagator survive, we should impose the scaling dimension of space and time as: 

[L] = a, [S] = a, (3.1) 

where a is the lattice spacing for space. From now, we shall extend the non-critical SFT based on GCDT 
without changing the scalings above. Since we think that the causality is the identity of CDT, this sort 
of extension is meaningful to get some deep understanding of what CDT is. 

First, we consider the strings with different charges, (+)-type and (— )-type. The creation and anni- 
hilation operators for (+)-type string, \I/+(L) and ^ + (L), and for (— )-type string, ^_(L) and ^-(L), 
are defined as the following vacuum conditions, respectively: 

W\.(L) = ^+(L)|0) = (0|^1(L) = V_(L)|0> = 0. (3.2) 

We assume these operators obey the following commutation relations: 

[^ + (L),^ + {L')] = ty_{L)rf_{L')\ = 6(L - L'), (3.3) 

and the others are zero. Additionally, we assume the same scaling dimensions with those of GCDT: 

[4(L)] = a°, VP±{L)]=a-\ [G s ] = a~ 3 , (3.4) 

where G s is the string coupling constant as before. Under the conditions above, we can extend the 
Hamiltonian for GCDT applying the interaction for spin clusters introduced by Ishibashi and Kawai [6j. 
Here we call such an interaction the IK-type interaction. It is based on the so-called peeling procedure in 
a discrete random surface. For example, considering a randomly triangulated surface coupled with Ising 
spins with one boundary and furthermore assuming that the boundary triangles have homogeneous spins 
(all spins are up- type or down-type), one peels triangles along with the boundary as if one peels an apple. 
If one continues to peel off triangles over the boundary triangles and one encounters the triangle having a 
different spin, then one surrounds the triangles having different spins by the triangles having same spins 
with the boundary triangles. In short, the randomly triangulated surface is separated by domain walls. 
In this case, the SDE for their approach coincides with the loop equation for the chain-type two-matrix 
model describing the random geometry coupled with Ising spins. We emphasize here that the above 
closed strings are not seen as the spin boundary as in the case of IK but the equal-time hypersurfaces 
with different charges. If we apply the IK-type interaction, we can write down the extended Hamiltonian 
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for GCDT: 



H m = 



dLil>\{L)U Q {L,k)ip + {L) + G s dLi dL 2 ipl{L 1 )i>l(L 2 )iP + (L 1 + L 2 )(L 1 + L 2 ) 



bG s 



dL x dL 2 ^\{L x +L 2 )^_{L 2 )^ + (L x )L x 
Jo 

/>oo 

dLi / dL 2 iPl(L x + L 2 )ij+(L 2 )^ + (L X )L 2 L X 



dL8{L)iP + {L) 



(3.5) 



where a and 6 are dimension-less constant^ In the Hamiltonian above, for simplicity, we will restrict 



G. 



olG. h 




Propagator Branching interaction Merging interaction 



Tadpole 



IK-type interaction 



Figure 2: Terms in the extended Hamiltonian: The red string stands for the one having the (+)-type 
charge, and the blue for the (— )-type charge. Of course, terms whose charges are flipped exist in the 
Hamiltonian, but here we do not include the graphical expressions of those terms. 



the topology of geometries to that of a disk, which can be realized by the following Hamiltonian: 



H» = lim H„ 

a— >0 



Next we will derive the SDE in our extended model. The SDE corresponds to Wheeler-DeWitt's equa- 
tion for the wave function of the universe. To begin, we define a partition function and disk amplitudes: 



Z = lim (Ole-^IO) = 1, 



S-voo 



and 

The SDE for W±(L) is 



W ± (L)= lim (0|e- SH "4(L)|0). 

o— >oo 



d 



lim ^(0|e- s ^4(L)|0)=0. 



(3.7) 
(3.8) 

(3.9) 



Using the equation, H^\0) — 0, and the commutation relations (3.3), we can rewrite the SDE as: 
= -Ld 2 L W ± (L) + ALW± (L) - 5{L) + G S L / dL x lim 



SH°h\ 



S— >oo 



^i(L!)^(L-Li)|0) 



bG s L I dL x \mv {Q\e- SH "^{L + L x )^Ul + L X )\Q). 



S— >QO 



(3.10) 



3 In fact, it is possible to include the interactions, f dLip^_(L)'Ho(L, A)ip+(L) and its spin-flipped term. However, because 
of the ^-symmetry as to the spin reflection, such terms merely cause a constant shift of the string coupling constant, so 
that we have not included these terms in the Hamiltonian. 



G 



Here we introduce the factorization theorem: 



lim ( 

S— >oo 



-SHR m iA 



^(i!)^(i 2 )|0) = lim 



S-s-oo 



^(£i)|0) Jim (0|e- Sff -4(i 2 )|0) 



(3.11) 



Applying the above factorization theorem, the SDE (3.10) becomes 



poo 

= -LdlW±(L) + KLW±{L) - 8{L) + G S L / dL 1 W±{L 1 )W±{L - L x ) 

Jo 

/>oo 

+ bG s L / dL 1 W±(L + L 1 )W^{L 1 ). 
Jo 



(3.12) 



In fact, our system has Z 2 -symmetry w.r.t. a spin-reflection, so that we will focus on a Z 2 -symmetric 
solution of the SDE: 

W+(L) = W-(L) = W A (L). (3.13) 

Next, we implement the Laplace transformation, £[Wa(£)] = J °° dLe~ LZ W A (L) = W A (Z). Applying 
the expression, W A (L), and Laplace transforming (3.12) yields 



= d, 



(Z 2 - A)Wa(Z) - G s W A {Zf 



-l + bG s £ 



L / dLiWAiL + L^Wa^) 



(3.14) 



We notice that the last term includes a divergent part as Z — > oo. To regularize this divergence, it is 
good to symmetrize it w.r.t. the reflection, Z f-> — Z [5] [7]: 



/•OO POO 

/ / d J Lie- z(i+ - Ll) M/ A ( J L + J Li)e +zil W r A(Li) + (Z O -Z) = W A (Z)W A {- Z) . 

Jo Jo 



(3.15) 



Subtracting the SDE with the reflection (Z — » — Z) from the original SDE (3.12|, we get the finite SDE 



= 8 Z 



(Z 2 - A) (W A (Z) + W A {-Z)\ - G s (W A (Z) 2 + W A (-Z) 2 + bW A (Z)W A (-Z) 



(3.16) 



Integration of the SDE above over Z yields 

c=(Z 2 - A) (W A (Z) + W A {-Z)) - G s (w A (Z) 2 + W A (-Z) 2 + bW A {Z)W A (-Z)) , (3.17) 
where c is a constant. 

We will derive a perturbative solution for the SDE above around the weak coupling region, G s < 1, 
by expanding the loop amplitude W A (Z) and c like: 



W A (Z) = J2GsW n (Z), c = J2g?c 



n=0 



n=0 



As for W (Z), we find 



W (Z) = 



1 



(3.18) 



(3.19) 



Z + VA 

where we have chosen an overall constant for Wo(Z) to coincide with that of pure CDT pQ. As for W\(Z) 
and Wi(—Z), we find 



Wi(Z) + W^-Z) 



1 



(Z + VA) 3 (Z- VA) 3 



ciZ 4 + (2 - b - 2Aci)Z 2 + ciA 2 + 2A + bA 



(3.20) 



Assuming that the disk amplitude behaves as 1/Z in the large Z-region, we can determine that c\ = 
— (6+ 1)/2A. Furthermore, we can extract W\{Z) by considering that Wi(Z) is analytic in the region, 
Rc[Z] > 0. Thus, the perturbative solution is 



W A (Z) = 



Z 



1 -G 1 

/A S 4A 



Z + 3vA 

JzWW 



(Z + VA) 



-0(G 2 S ). 



(3.21) 



The solution with b — is equivalent to that of the pure GCDT (2.7). 



where 



3.2 Matrix Model Approach 

We start with the following matrix integral: 

1> + d<l>-e-£ v U+'*-\ (3.22) 

0_) = - 9 (0+ + 0-) + i(0 2 + + 4>-) - §(</>+ + 4>l) + x<t> + 4>-. (3.23) 

In the integral above, <j>±, g, g s and x are N x N hermitian matrices, the 'tHooft coupling constant, the 
string coupling constant and the coupling constant characterizing the interaction, respectively. Then, we 
expand the fields and coupling constants w.r.t. the lattice spacing a as follows: 

(p+ = I - a(A + B) +0{a 2 ), 0_ = J — a(A - B) + 0(a 2 ), (3.24) 

and 

g s = a 3 G s , g= ^(l- ^a 2 (A - 2A) + 0(a 4 )^ , x = Xa 2 , (3.25) 

where A and B are N x N hermitian matrices, and I is the unit matrix, and G s , <3?, A and X are the 
corresponding renormalized values. Thus, our model can be seen as the one that the cut-length shrinks 
to zero (<7s —> 0), and the strength of the interaction falls off (x — > 0). The causality induces the scaling, 
g s — > 0, and in addition, by taking the limit, x — > 0, we can get our model as the weakly interacting 
model. Substituting the fine-tuned values, we can write down the partition function of the matrix model 
in the continuum limit: 



dAdB exp 



~tr( ^ A 3 + AB 2 -AA 



(3.26) 



An interesting thing is that in the matrix model having this type of potential, the Gaussian integral over 
B can be performed by introducing the eigenvalues Aj's for the matrix A [S]: 

Zoc yn^ 4 A 2 (A)r[(^ + A,r 1/2 e-* y , (3.27) 

i i,j 

where 

N N, v 
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and A(A) denotes the Vandermonde determinant, A(A) = Ili<j(^i ~ ^i)- ^ n the large- N limit, the saddle 
point equation becomes 

A J 



G, 



where V'(Xi) = A 2 — A. Here we define the resolvent fo r A as W A (Z) = -^tr(Z— A)" 1 , and the distribution 
of eigenvalues as p(X) = jj J2 t (5(A — Xi). Multiplying (3.29 1 by 1/(Z— Xi) and summing over i, we obtain 
the loop equation in the large- TV limit: 

V'(Z)W A (Z) + V'(-Z)W A {-Z) - G S (W A (Z) 2 + W A {Z)W A (-Z) + W A (-Z) 2 ) + G s n(Z) = 0, (3.30) 



where 



G s n(Z) = / dXp(X) 



V'(X)-V'(Z) V'(X) - V'(-Z) 



Z - X 



dXp(X)X. 



Z + X 



(3.31) 



In the calculation above, we used the two identities: 

1 



2 



Z-^ 7, - 



. . Z — Xi Xi — Xj 



W A {Z) 2 



~Wi(Z), 



(3.32) 



and 



-y— 



TV 2 



1 



1 



\$ -\- \j V Z — X[ Z -\- \^ 



= -W A (Z)W A {-Z). 



(3.33) 



Putting explicit form of the potential into the loop equation (3.30), we find 
(Z 2 - A) (W A (Z) + W A {-Z)\ - G s (W A (Z) 2 + W A (Z)W A (-Z) + W A (-Z) : 



dXp(X)X. (3.34) 



Remembering the SDE derived in the non-critical SFT approach (3.17), we can find a great similarity 
between the two. Namely, if we set b = 1 in the SDE, then the two equations are exactly same. Thus, 
this matrix model in the continuum limit can reproduce our GCDT with extended interactions in b = 1. 

We can extend the matrix model in the continuum limit above to the general 0(n) vector model [5] 
such that: 

Z= [ dAdB 1 ---dB n e-^ u{A B ^ A) , 



where 



U(A,Bx, ■ ■ ■ ,B n )=A(Bf 



Bl) 



-A 6 - AA. 



(3.35) 



(3.36) 



and A, B\, . . . , B n are N x N hermitian matrices. One can find that the previous matrix model in the 
continuum limit is 0(1) vector model. Again, we can integrate out all Bi's, and a consequence is 



r N 

i=l i,j 



(3.37) 
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where A^s are eigenvalues of A, and V = £^ — J^j(§Af — AAi). A saddle-point equation becomes 



2 

TV 



(3.38) 



In the similar manner as O(l) vector model, we get the loop equation for the resolvent W A (Z): 

(Z 2 - A)(W A (Z) + W A {-Z)) - G S (W A (Z) 2 + nW A (Z)W A (-Z) + W A {-Z) 2 ) = 2 J d\p(X)X. (3.39) 

Thus, the loop equation of this 0(n) vector model coincides with the SDE labeled by a free parameter b 
(3.171 only if we identify n with b. 

3.3 Inclusive Process 

In the above, we derived the differential equation for disk amplitude in our extended model, and solved 
it by perturbative expansions. As a confirmation, we shall reproduce the same differential equation for 
disk amplitude using the so-called inclusive process ([5], [5] and [H]). In the inclusive process, putting 
caps (disk amplitudes) on one of two loops (universes) at branch points we can focus on the amplitude 
with one loop, which has its origin in the so-called transfer matrix formalism [13] . If we focus on the case 
that initial and final strings have the same charges, then the inclusive SFT Hamiltonian can be written 
as follows: 



Hin — 



dL^ + (L)H (L, k)ip+{L) + 2G S 



dLi 



dLiW+iLMliLity+iL! + L 2 )(Lt + L 2 ) 



bG s 



dh x 



dL 2 ^UL 1 + L 2 )W_(L 2 )^ + (Li)Li + 



Cylinder amplitudes are denned as: 



pOO pOO 

G ++ (L U L 2 ) = / dSG ++ {L 1: L 2 -S)= / (0|^ + (L 2 ) e - s ^«Vl(ii)|0), 
Jo Jo 

poo poo 

G__(Li,La)=/ dSG—{LiM;S)= / (0|^-( J L 2 ) e - 5 ^«^L( J L 1 )|0). 
Jo Jo 

By a differentiation of G++(Li,L 2 ) w.r.t. time S, we find 

d s G ++ {L x ,L 2 ;S) = -(0|V;+(£ 2 K SHj "[#/iV,^( L i 



(3.40) 

(3.41) 
(3.42) 

(3.43) 



where we have used Hin\0) — 0. With the calculation similar to (3.10), the equation above can be 
rewritten as: 



(3.44) 



dsG+^LuL*; S) = L x {d l Li - A)G++(L 1 ,L 2 ; S) - 2G s L l / dLW+{L x - L)G++(L, L 2 - S) 



bG s Li J dLW-(L)G ++ (L 1 +L,L 2 ;S). 



Limiting the length L 2 to and integrating over S in (3.44), we have 



= Li{d 2 Ll - A)W+(L 1 ) - 2G S L 1 J dLW+{L x - L)W+(L) - 6G s Li / dLW-(L)W + (L x + L), (3.45) 
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where W+(L) = J °° dSG ++ (L, 0; S), and we have used the fact that G++(Li, 0; oo) = G++{L\, 0; 0) = 0. 
Our system has the Z2-symmetry as to the spin reflection, so that we focus on a Z2-invariant solution, 
Wa(L) = W±(L), as in (3.13). Then, implementing the Laplace transformation of (3.45) yields 



= d x 



-{X 1 - A) + 2G S W A {X)\ W A (X)+bG s / dLxe 



-XL! 



dLW A (L)W A {Li +L) 



(3.46) 



where W A (X) = / °° dLe LX W A (L). Again, the last term includes divergent part as X —> oo. Thus, we 
;he diverg 

(X 2 -A) (w A (X) + W A (-X) 



need to remove the divergence by the symmetrization as in (|3.15) 

W A (X) 2 



= d x 



2G, 



~W A (X)W A (-X) 



W A {-X) 



(3.47) 



After the proper shifts of the string coupling constant G s and the free parameter 6, the equation (3.47) 
coincides with the SDE for the disk amplitude (3.17) as expected. 



Next, we start with a discrete model, and then reconstruct our model as its continuum limit. To carry 
it out, based on the transfer matrix formalism 13J we derive the differential equation for disk amplitude, 
which turns out to be equivalent to (3.17) in the continuum limit. Here the transfer matrix is the one- 



time-step propagator having the length-Zi initial loop with (±)-charge and the length-^ final loop with 
(±)-charge denoted by G±±(/i, h] !)■ First, we derive the non-interacting propagator which is one of parts 
in the transfer matrix, G±(ii, hi 1)- This can be easily calculated introducing the generating function of 
it, i.e. &±±{x±, y±; 1) = l2 x±y±G±±(li,l2\ 1). Namely, we prepare four types of triangles weighted 
by gx± and gy± (Fig. [3]), and only from the combinatorics we can find the generating function of the 
one-time-step propagator [T|: 




gx + 




gy+ 




gx- 




gy- 



Figure 3: Four types of triangles weighted by gx± and gy±: gx± and gy± are the weights for the 
triangles which are the elements of an initial loop with (±)-charge and of a final loop with (±)-charge in 
the generating function, respectively. 



G r ±± (x ± ,y ± ;l) = (3.48) 

- gx±)(l - gx± -gy±) 

In the calculation above, we marked a point on one of initial links following pQ. For the later discussion, 
we give the specific form of weights: 

g= 1 -e-^ A , x ± =e~ aX ±, y = e~ aY \ (3.49) 

where a is the lattice spacing, and A, X± and Y± are the renormalized coupling constants. Then, wc 
introduce the transfer matrix combing the disk amplitude with the length-/ initial loop with (±)-charge 
w±(l) and the non-interacting propagator Gj_±(/i, Z2; 1) as follows: 

00 

G ±± {h,l 2 \ 1) = Gl£h,h\ 1) + 2g s J2 hw ± (h - l)Gi ± {l, l 2] 1) + bg s £ hw^G^l + h,l 2 ; 1), 

1=1 1=1 

(3.50) 
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Figure 4: Terms in the transfer matrix, G++(Zi, I2', 1): dots in the arguments are replaced by some suitable 
variables on a case-by-case basis, and dots on loops in pictures are marked points on links. 



where g s is the bare string coupling constant and & is a free parameter (Fig. |4J). A natural property of 
the propagator is the decomposition law: 



G ±± (h,l 2 ; s) = J2 G ±± (h,l; 1)G±±(Z, h\ a - 1). 



1=1 



After the Laplace transformation, the equation (3.51) becomes 

/ dz± 



G±±{x ± ,y ± ;s) = j> G±±(x±, z ± 1 ; l)G±±(z±, y±; s - 1). 



(3.51) 



(3.52) 



Substituting (3.481 into the equation (3.52 1, one finds 

dz± 



G±±(x ± ,y ± ;s) = j> ^g- ]T x^z^G^h, I; 1)G ±± (Z, l 2 ; s - 1) 



J 27TIZ± 



il,/3,i=l 

G± ± (x±, 1) + 2g s x±d x± ( u> ± (x±)G± ± (x±, z± x ; 1) 



j,=i ;'=i 



E E ^^±^± ( 4^t(^)G4 ± (/' + Ix,^ 1 ; i; 



G±±[z±,y±;s-l), (3.53) 



where w)±(x±) = J^i x l ±w± (I), and G T ±± (l' + h,z ± 1 ;l) = E,i^± Gl ±± (l' + h,l;l)- Here we introduce l C ut 
to regularize the divergent summation over I' in (3.53): 



G±±(x±,y±;s) 



1 + 2g s x±(d x± w±(x±) + w±(x±)d±) 



gx± 



-G 



±± 



fJ 



OO Z cut 

ii=i ('=1 



2iriu± 



G±± 



<fa± -(i'+ii) ff M ± 



1 - gu ± 



y±\s- 1 
a 



1 - 0U± 



, y± : s - 1 



(3.54) 



where u± = e~ aU± . 

In the following, we focus on Z2-symmetric solutions, i.e. G\(x,y; s) = G±±(x±,y±; s) and w\(x) = 
w±(x±). Under the scalings, i.e., S = as, L\ = al\ and L2 = (1I2, one finds the following renormalized 
functions: 



G A (X,Y;S) = lim aG x (x,y;s), W A (X) = lim aw x (x), W\(L) = lim w\(l). (3.55) 

a— *-0 a— »0 a— J"O0 
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From (3.49), (3.53), (3.54) and (3.55), one finds 



dsG A (X,Y;S) = -d x 



{X 2 - A)G A (X, Y; S) + 2G S W A (X)G A (X, Y; S) 



lim bG s / dLx / dL'e W A (L')e 



G A (L' + L 1: Y;S) 



, (3.56) 



where L cnt = al cut , and we have used G A (V +Lx,Y;S) = dUe( L ' +Ll ^ u G A (U, Y; S). Implementing 
the inverse Laplace transformation w.r.t. Y and integrating over S in (3.56), one finds 



= -d 



'x 



(X 2 - A)G A (X, L) + 2G S W A {X)G A {X, L) 



bG s I dLx / dL'e L ' x W A (L')e 



-Is -(L'+LJX 



G A (L' + L U L) 



(3.57) 



where G A (X,L) = f n dSG A (X, L; S), and G A (L' + L\,L) = f™^ dY e LY G A (L' + L U Y) . Limiting the 
length L to in (3.57), one finds 



o = d x 



(X 2 - A)W A {X) + 2G S W A {X) 2 + bG s I dLx I dL'e L ' x W A (L r )e 



W A {L' + Lx) 



(3.58) 



where W A {X) = G A ( X 0) a nd W A {L' + L x ) = G A {L' + L x , 0). As for the last term in ( |3.58[ ), we use the 
same procedure as in (3.15) and ( |3.47[ ): 



/•oo />oo 

/ dLx dL'e L ' x W A {L')e^ L ' +L ^ x W A {L' + Lx) + {X ^ -X) = W A {-X)W A {X). (3.59) 
"'O Jo 



Therefore, we obtain the finite differential equation for disk amplitude: 



= d x 



(X 2 - A) (W A (X) + W A {-X)) + 2G S ^a(^) 2 + \w A {X)W A {-X) + T^a(-A^) 2 



(3.60) 



Finally, after the proper shifts of the string coupling constant G s and the free parameter b, one finds that 



(3.60) is equivalent to (3.17) as expected 



4 Discussions 

We have shown the equivalence between the two different field theories at the level of differential equations, 
the Schwinger-Dyson's equation in the non-critical SFT and the loop equation of the matrix model in the 
continuum limit. We hope that our model is a first step toward matter-coupled systems based on CDT. 
In the following, we will examine the model constructed in this paper from different point of view. 

To begin with, we will discuss our model in terms of the SFT approach. Although we have used the 
IK-type interaction to construct the extended SFT based on GCDT, we do not understand whether or 
not our model is on the critical point of the Ising model characterized by the Curie temperature. In the 
following, we will explain two complications around this problem. First, at the critical point of Ising 
spins the spin configuration must be random. In other words, the spins are supposed to fluctuate all 
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length scales between the lattice spacing and the correlation length. Contrary to that, in our model 
the homogeneous spin (charge) configurations survive as the propagators. Second, the definition of time 
induced by our Hamiltonian ( |3.5[ ) is different from the would-be GCDT coupled with Ising spins. Namely, 
we consider the closed strings in our model as not spin-cluster boundaries but spacial boundaries, so that 
we pursue the time flow of spatial boundaries. Thus, our time is nothing but the proper time. This 
proper time is crucially different from the time defined via the spin-cluster boundary [TJ] [TS]. If we 
consider our time as the one defined via the spin-cluster boundary, which is equivalent to treating our 
model as the GCDT coupled with Ising spins, then the scaling dimension of time may be different from 



the lattice spacing a according to [14]. This contradicts our first setup (3.1 1. Anyhow, the free parameter 
b might be the key to know what our model is. 

In addition, it is possible to extend our non-critical SFT to the multi- "colored" system: 

n poo n poo poo 

H £ )= Y, dL4(L)H (L,A)ip i (L) + G s y2 / dL x dL 2 i/jl{L 1 )ipl(L 2 )ip l (L 1 + L 2 ){L 1 + L 2 ) 
i=l J o i= i J o J o 

n n pOO pOO 

+ G *y2J2 b v dLi di2^(ii+i 2 )^(i 2 )^(ii)ii 

i=l jfr J o J o 

n pOO pOO 

l=1 Jo Jo 

n poo 

-J2 dL6(L)ipi(L). (4.1) 

We can derive the free parameter b in our model from the multi- "colored" system above under the 
treatment, W\{L) = ■ ■■ W n (L) = W\(L), bij = for j = i and b^ = 1 for j ^ i. 

Next, we will closely look at our matrix model. Considering the direct product of the two copies of the 
potential, each of which yields the pure GCDT, and introducing the linear combinations of the matrices 
as $ + = A + B and $_ = A - B, we find 

i Q$ 3 + - A$ + + V - A$_) = i- Qa 3 + AB 2 - AA) , (4.2) 

where G s — 2G S . This is the potential of our 0(1) vector model in the continuum limit. Then, diago- 
nalizing the matrix A as A — diag(Ai, • • • , Xn) and integrating out the matrix B, we get the effective 
theory for the eigenvalues of A with the potential, 



-— (terms induced by the integration over B). (4.3) 



terms appeared in the pure GCDT 



The important point here is that our model is slightly different from the pure GCDT matrix model 
because integrating out the matrix B an extra correction is added to terms appeared in the pure GCDT. 
From the matrix A's point of view, the matrix B can be seen like some external field. The strength of 
such an external field can be bigger by inserting the integrated-out matrices, which leads to the 0(n) 
vector model in the continuum limit. 
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